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ABSTRACT 

A comprehensive CFD-based modeling approach is 
applied to several automotive catalytic converters, with 
the objective of predicting thermal behavior under 
steady-state, high-load conditions. Specialized computa- 
tional models are used to account for effects of heat and 
mass transfer in the monolith, oxidation reactions, heat 
generation, conjugate heat transfer in the various con- 
verter materials, and radiation heat transfer. These vari- 
ous physical considerations are assembled in a 
comprehensive CFD model, which is solved using state- 
of-the-art computational techniques. Detailed tempera- 
ture measurements, taken in engine-cell experiments at 
Ohio State University's Center for Automotive Research, 
are used to validate the CFD models. Excellent agree- 
ment seen between measured and computed tempera- 
tures, both inside the converter assembly and on the 
outer shell. Trends in the predictions and experiments 
are discussed, with the aim of understanding the physical 
mechanisms of heat transfer in monolithic catalytic con- 
verters. An increased understanding of these heat trans- 
fer mechanisms can be used to improve converter 
durability at high temperatures. 

INTRODUCTION 

Temperature behavior is an important consideration in 
exhaust system design. Numerical methods, including 
computational fluid dynamics (CFD), are often employed 
to predict temperatures of exhaust systems, including 
catalytic converters. Temperature predictions are used to 
predict mounting system durability, to assure the viability 
of converter location relative to heat-sensitive compo- 
nents, and as input to structural finite-element analysis. 
In catalytic converter flows, numerical techniques are 
made more complex by the presence of chemical reac- 
tion and heat generation in the converter monolith. In 
addition to chemical reaction, catalytic converters include 
such features as porous media flow, conjugate heat 
transfer between materials of various properties, and 
thermal radiation heat transfer at high temperatures. 



Specialized computational methods are required to deal 
with these effects in a numerical simulation. 

Thermal studies of catalytic converters have been per- 
formed in the past, often with the objective of predicting 
emissions and thermal performance in transient situa- 
tions. Early work by Young and Findlayson (1) and Sher- 
ony and Solbrig (2), among others, treated the problem of 
heat and mass transfer in monolith channels using ana- 
lytical methods. Voltz et al. (3) first proposed oxidation 
rate equations for carbon monoxide and propylene on 
platinum-washcoated monolith channel walls. These 
equations are used in a variety of converter light-off stud- 
ies. Using these reaction rate expressions, Heck et al. 
(4), among others, investigated heat & mass transfer with 
reaction kinetics in a single monolith channel. Oh and 
Cavendish (5) extended the method to investigate cata- 
lyst overheating problems under transient driving condi- 
tions. 1-D methods based on the same reaction rates 
and heat & mass transfer relationships are used exten- 
sively in the literature; see for example Siemund et al. (6) 
and Schweich (7). 

The results obtained in single-channel studies were 
extended to two spatial dimensions, including the effects 
of a non-adiabatic wall boundary, by Chen et al. (8). The 
computed results were used to predict severe thermal 
gradients in the monolith, which have been shown to lead 
to monolith cracking problems. In a similar study, the 
method was extended from platinum-only catalysts to 
platinum/rhodium "three-way" catalysts by Chen and Cole 
(9). The multi-dimensional nature of these studies make 
them useful for predicting temperature fields, which can 
have a significant impact on mat erosion, shell deforma- 
tion, and monolith cracking. 

In the above-listed cases, constitutive equations are 
applied only to the monolith region, and not to the 
exhaust gas flow upstream and downstream of the cata- 
lyst. Velocity distribution across the monolith cross- 
section has been either an assumed function or assumed 
to be uniform at the monolith inlet. The relationship 
between monolith flow distribution and total pressure loss 
has been studied in some detail however, experimentally 
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by Wendland et al. (1 0, 11 ), and computationally by Ben- 
jamin et al. (12), Weltens et al. (13) and Lai et al. (14). 
These studies investigate the effects of inlet cone angle 
and substrate geometry on such performance parame- 
ters as monolith flow distribution and backpressure. 

Generally speaking, the flow-related computational work 
in the open literature ignores heat generation effects of 
oxidation reactions. Conversely, computations which aim 
to predict reaction kinetics and exothermic behavior in 
converter monoliths often neglect effects of flow outside 
the monolith itself. In order to fully understand the conse- 
quences of converter design, a model must be formu- 
lated to include both sides of the problem; both reaction- 
based heat transfer and fluid flow. 

The current work incorporates a fully three-dimensional 
treatment of all the relevant flow and heat transfer phe- 
nomena in catalytic converters, including; 

• heat & mass transfer between the monolith solid & 
exhaust gas; 

• oxidation reactions of carbon monoxide, propylene, 
and methane, and associated heat generation; 

• 3-D Navier-Stokes equation treatment of flow 
upstream and downstream of the monolith; 

• laminar anisotropic porous media treatment of flow in 
the monolith; 

• 3-D, fully conjugate heat transfer in the monolith, 
insulating mat, and converter shell; 

This comprehensive modeling methodology is applied to 
three different converter designs. Each design is simu- 
lated at two different high-temperature, steady-state 
engine conditions, such as would be used in converter 
durability tests. The results are compared to correspond- 
ing engine cell experiments in order to validate the com- 
putational method. Of special interest here are 
temperature fields in the catalytic converter. Predicted 
and measured pressure loss and conversion perfor- 
mance data are also presented. 

PHYSICAL MODELS 

POROUS MEDIA FLOW -The governing equations of 
flow and heat transfer in a catalytic converter can be 
divided into two general categories: those equations 
which apply in the monolith brick itself, and those which 
apply in the regions outside the monolith. In the monolith, 
flow is assumed to be laminar and unidirectional along 
the converter axis (defined here as the z-axis in Carte- 
sian coordinates). Neglecting entry effects, pressure 
loss in a monolith channel is given by Shah and London 
(15) as 



AP_ 2^(/Re) 1 



^ 7 = (/Re)(4z + ) 



which can be rearranged as 
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The friction factor f times channel Reynolds number Re is 
a constant quantity which changes only with cross-sec- 
tional shape. For a square-channel cross section, (fRe) = 
14.227 (ref. 15). 

MONOLITH ENERGY CONSERVATION - The energy 
conservation equation for the fluid is related to the energy 
conservation equation for the solid using the local source 
term 



S x =ha s {T 5 -T f ) 

such that in the fluid at steady-state; 
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Note that equation 4 does not consider heat diffusion 
across individual fluid channels in the x and y directions. 
This simplification is often made (8,9) because the diffu- 
sion is dominated by the surface interchange term on the 
right hand side of equation [4]. Put another way, the 
heat flux between solid and fluid described by equation 
[3] is far greater than the heat flux across a fluid channel 
diameter: the latter is therefore neglected. 

In the monolith solid, the heat conduction equation is 



y 

i=3 



(5) 



Note that the solid conductivity X is anisotropic due to the 
monolith structure of the solid. 

CHEMICAL REACTION AND HEAT RELEASE - The 
final term in equation 5 represents the heat of oxidation 
reaction for carbon monoxide (CO), propylene (C 3 H 6 ), 
and methane (CH 4 ). Carbon monoxide is well known to 
be a prime contributor to heat generation in the monolith. 
The combination of propylene and methane is meant to 
simulate the presence of unburned hydrocarbons. These 
two species have been used in the past to represent 
unburned hydrocarbons (ref. 5), recognizing that many 
hydrocarbon species are typically present in real exhaust 
gas. The oxidation of these species is assumed to take 
place on the surface of the monolith channels, once the 
reacting species have been transferred from the gas 
phase to the solid phase. Reaction rates Rj are functions 
of local temperature and surface concentrations 



(6) 
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where 



CFD METHODOLOGY 



G = T(l + K,c m + K 2 c C3HS f (l + K 3 c 2 C0 c 2 C3H J\ + K,c% ) 

k 1 =k co = (6.699x1 0 13 )exp(-12556/T) 

k2=k C 3H6= C 329x10 15 )exp(-14556/T) 

k 3 =k CH4 = (7.326x10 10 )exp(-19000/T) 

K n = 65.5exp(961/T) (6A) 

K 2 = 2080exp(361/T) 

K 3 =3.98exp(11611/T) 

K 4 = 479000exp(-3733/T) 

MONOLITH HEAT AND MASS TRANSFER - Mass 
transport equations, analogous to the thermal balance of 
equation 4, are used to represent species diffusion from 
the gas to the solid in the channels; 

V z ^- = m mi a s (c fl -c si ) (7) 
dz 

In equations 4 and 7, local heat and mass transfer coeffi- 
cients (h in equation 3 and m m in equation 7) can be esti- 
mated using one of several existing relationships. Shah 
and London (15) give the following relationship for lami- 
nar flow in square cross-section channels; 

Nu =^ = 2.97 (8) 

which can be used to evaluate the channel heat transfer 
coefficient h. Equation 8 can be extended to the corre- 
sponding mass transfer relationship; 

5/1 = ^^ = 2.97 (9) 
D 

^ mi 

which can be used to find mass-transfer coefficient m mi . 

It should be noted that, although the reaction rate expres- 
sions are more complex, it will be the Sherwood number 
which limits the reactions in the cases performed here. 
At high temperatures, the reaction is said to be mass- 
transfer limited, and will proceed almost instantaneously 
once species transfer has been made from the fluid to the 
solid phase. The constant values used here are chosen 
for their theoretical soundness and ease of implementa- 
tion. Though more detailed expressions for Nu and Sh 
are the subject of much investigation (16,17), the deter- 
mination of such expressions are beyond the scope of 
this report. 

Outside the monolith, flow is governed by the Navier- 
Stokes equations in three dimensions. The heat conduc- 
tion equation is used to represent the heat conduction in 
the shell and insulating mat. 



SUMMARY OF CASES - Three catalytic converter 
designs are chosen for the study, as listed in Table 1. In 
case 1, the catalytic converter monolith is not wash- 
coated, and therefore does not include significant chemi- 
cal reaction. (In the experiments, an active, washcoated 
converter is inserted upstream of the test converter in 
case 1 in order to ensure that no reactants are present in 
the test piece). Case 1 is meant to evaluate the heat- 
transfer model without any reaction or exothermic behav- 
ior. Case 2 uses the same converter geometry as case 
1, but in case 2 the monolith has been treated with a 
platinum-impregnated washcoat. Case 2 will therefore 
include significant chemical reaction and heat genera- 
tion. In case 3 the converter design is significantly differ- 
ent. There are now two separate monoliths of oval cross- 
section, and there is no insulating mat in the cones or 
between the two monoliths. The washcoat in case 3 is 
equivalent to the washcoat in Case 2. 

Each case listed above is simulated and tested at two 
engine conditions in order to thoroughly validate the mod- 
eling methodology across a range of temperatures. The 
two engine conditions are referred to as M for medium, H 
for high. So for example, the case 2 converter using the 
high engine condition will be referred to as case 2H. Vari- 
ous cases performed are summarized in Table 1. 
Sketches showing the internal layout of the converters 
are shown in Figure 1 . 

The converter materials are typical of current technology 
in use in the automotive industry, including vermiculite- 
based expanding mat mount material, 409 stainless steel 
shell, and ceramic 400cell/in 2 substrate. Because the 
reaction rate expressions given by equation [6] are for a 
platinum-only catalyst coating, and because this is a vali- 
dation exercise, the coating on the active substrates is 
platinum-only, as opposed to a more typical platinum/ 
rhodium or palladium formulation. (Expressions for Pt- 
Pd-Rh combinations to be used in the place of eqn. [6] 
are scant in the literature, and development of such 
expressions are beyond the scope of this report.) 



Table 1 : Summary of Cases 





cross- 


number of 




engine 


case 


section 


monoliths 


washcoat 


condition 


1M 


round 


r 


none 


medium 


1H 


round 


V 


none 


high 


2M 


round 


V 


35 grit 3 Pt 


medium 


2H 


round 


V 


35 gm 3 Pt \ 


high 


3M 


oval 


2 


35 gm 3 Pt 


medium 


3H 


oval 


2 


35 gm 3 Pt 


high 



* two monoliths adjacent to each otherp.e. butted-up) 
will be treated as a single monolith 
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COMPUTATIONAL MODELS -A computational grid is 
generated for each of the six cases to be studied using 
ICEM-CFD's Hexa grid generation software. Each grid is 
composed of approximately 310,000 hexahedral cells, 
representing a quarter-symmetry geometry of each con- 
verter. The mat and converter shell are both represented 
by 3-D hexahedral cells. The grids are shown in Figure 2. 
Note that Figure 2 shows only the part of the grid com- 
posing the converters: the grid extends upstream , out of 
the range of Figure 2, to represent the inlet pipe section. 

The equations of flow are discretized using the self- 
filtered central-differencing scheme (SFCD). This 
second-order scheme offered the best combination of 
robustness and accuracy for this class of cases. The 
standard k-e turbulence model is used to account for tur- 
bulent fluctuations. In the near-wall region, the low Rey- 
nolds number k-e model of Norris and Reynolds (18) is 
applied. This "two-layer" turbulence modeling approach 
is often more accurate in representing near-wall turbulent 
flow behavior and has been recommended for use in this 
class of flows (12,19). 

The CFD practices listed briefly here are based on rec- 
ommendations of the literature, experience in CFD mod- 
eling, and a rigid interpretation of CFD guidelines. Every 
attempt was made to avoid common deficiencies in the 
open literature such as coarse or skewed grids, incom- 
plete computational domains, diffusive first-order discreti- 
zation schemes, and unsuitable turbulence modeling 
practices. The simulations here attempt to be of the 
highest quality in all aspects of general CFD practice 
before consideration of the specialized models for mass 
transfer, reaction, etc. 

User subroutines are required to perform the following 
tasks in the CFD simulations; 

• Solve for fluid and solid species concentrations in the 
monolith. At each cell location in the monolith, user 



subroutines equate the mass transfer rate (from fluid 
to solid phase) with the reaction rate in the solid 
phase (i.e., the elimination of reactant mass in the 
solid phase); 

• Generate heat in the solid phase, and transfer heat 
between solid and fluid phases in the monolith; 

• Calculate a local resistance to flow using the local 
temperature-dependent laminar viscosity and equa- 
tion [2]. 

The Star-CD user subroutines to carry out these opera- 
tions have been developed by Clarkson (20, 21). Theses 
subroutines, modified slightly where necessary are used 
in this study. 



internally 
insulated cones 
(inlet and outlet) 




' msuiavnamal ////// J 



flow 



monolith 



/ / / insulating mat /" 



a) cases 1 & 2 



flow 



Vy/ms. mal'/A 


\/ t ms. maV/n 


front monolith 




rear monolith 








Wins. ma\Y/\ 


V.ms.maV//\ 



no cone 
insulation " 



b) case 3 

Figure 1 . Catalytic Converter Geometry Cross-Sectional 
Views 




Figure 2. CFD Computational Grids 
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BOUNDARY CONDITIONS - For this type of validation 
study, it is imperative that the conditions assigned to the 
computations match the conditions seen in the experi- 
ments. Boundary conditions for inlet mass flow, tempera- 
ture, and species concentrations are therefore assigned 
to match experimental measurements as exactly as pos- 
sible. This boundary-condition matching leads to CFD 
boundary conditions which vary slightly from case to 
case, since the engine-cell experiments vary slightly from 
case to case. A summary of boundary conditions is 
given in Table 2. 

Table 2: CFD Boundary Conditions 



case 



inlet 


1M 


1H 


2M 


2H 


3M 


3H 


velocity (m/s) 


65.1 


82.8 


64.6 


116.8 


64 


107.2 


temperature (C) 


768 


803 


738 


829 


727 


824 


inlet CO (vol %) 






0.94 


0.93 


1.01 


0.31 


inlet total HC (ppm) 






677 


394 


942 


392 


inlet 02 (vol %) 






0.33 


0.18 


0.39 


0.21 


turbulence 


5% turbulence intensity Q 5mm length scale 


ambient wall 














film coeff (W/m2) 


20 


20 


20 


20 


20 


20 


temperature(C) 


303 


303 


303 


303 


303 


303 


shell emissivity 


0.3 


0.3 


0.3 


0.3 


0.3 


0.3 



inlet conditions based on measured data (turbulence is estimated) 



- constant inlet velocity and temperature profiles are 

allowed to develop over long computational inlet length 

- ambient conditions estimated based on straight-pipe temperature 

measurements 

MATERIAL PROPERTIES - A variety of approaches are 
required to calculate material properties. The gas proper- 
ties are dependent on temperature but not gas composi- 
tion, an assumption acceptable when the reacting 
species compose only a small fraction of the overall flow. 
This is referred to as a "passive - contaminant" 
approach. A summary of material property calculations 
and assumptions are listed in Table 3. 



Table 3 : Material Properties 



monouxn 




material 


— : 

cordierite ceramic 


cell density 


400cpst/6.5 thickness 


coating 


Pt-only 


gas 




rnai&nai 




density 


ideal gas law 


viscosity 


Sutheriands formula 


thermal conductivity 


polynomial function 


mass diffusivity 


reference (22), pg. 505 


molecular weight 


30.3 




(gas properties are not 




inTiuencoo oy reaction j 


shell 




material 


409 stainless steel 


thermal conductivity 


constant 27W/m A 2 


mat 




material 


intumescent expanding mat 


thermal conductivity 


polynomial fit (based on 




manufacturers data) 



Star-CD version 3.0c is the CFD software used to assem- 
ble and solve the problem. The solutions require approx- 
imately 800 iterations to converge. Approximate solution 
time is 36 hours on a Sun Ultra2. Table 4 lists several per- 
tinent computational parameters. 



Table 4: CFD Solution Parameters 



grid size 


310,000 cells 


grid type 


all hexahedral, structured-type grid 


discretization 


self-filtered central differencing 


solution algorithm 


SIMPLE 


convergence 


normalized residuals<1e-3 
(norm. residuals decrease 3 orders 
of magnitude for k, e equations) 


approx. run time 


36 hrs/case 
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EXPERIMENTAL METHOD 

Validation experiments were performed in an engine cell 
at Ohio State University's Center for Automotive 
Research. A sketch of the experimental system is shown 
in Figure 3. Note that in cases 1M and 1 H, where the test 
piece is not washcoated with catalyst material, an active, 
washcoated catalytic converter is placed in the exhaust 
gas stream upstream of the test piece. This is done to 
clean the exhaust gas of any species which might other- 
wise react in the test piece. In the remaining cases, the 
non-instrumented converter is replaced by a straight 
pipe. Experimental data was taken after the system had 
reached steady-state for a given engine condition. The 
following variables were recorded; 

• gasoline flow rate; 

• air/fuel ratio; 

• species concentrations of CO,02,C0 2 ,NO x , and 
hydrocarbons, upstream and downstream of the cat- 
alyst; 

• temperatures on the converter shell; 

• temperatures in the converter bricks; 

• pressure upstream and downstream of the catalyst; 



• gas temperature upstream and downstream of the 
catalyst. 

A diagram showing the system layout is shown in Figure 
3. The silencing element is inserted upstream of the con- 
verter to suppress severe pressure fluctuations which 
can affect the catalyst flow & heat transfer behavior. Tem- 
peratures were measured on the converter shells and 
inside the monoliths using K-type thermocouples. The 
locations for temperature measurements are shown in 
Figure 4 for the round converters (cases 1 and 2) and in 
Figure 5 for the oval converter (case 3). 

Some notation which will be used to describe the results 
is shown in Figures 4 and 5; 

• z/L - measured axially form the beginning of the con- 
verter shell, such that z/L = 0 at the beginning of the 
inlet cone and z/L = 1 at the end of the outlet cone; 

• z'/L - measured axially from the beginning of the 
front substrate, such that z7L=0 at the monolith front 
face and z7L=1 at the monolith rear face; 

• r/R mw - measured radially from the center of the con- 



verter such that r/R n 



0 at the monolith center and 



r/R mw =1 at the monolith outer radius. 




T 



I 



1 



V6 engine 




exhaust gas flow 



catalytic converter 
inserted here in cases IM. IH 



Emissions Measurement: 

CO. 0 3 : Horiba AIA-722 NDIR 

HC: Horiba FlA-720 



pressure 
measurement 



<\talytic converter^ 
test piece f 



exhaust 
gas out 



estimated uncertainty 




mass flow rate (velocity) 


+/- 


8% 


emissions species 


+/- 


1.5% 


temperature 


+/- 


0.75% (7-9 C) 



Figure 3. Experimental Exhaust System 
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zA=1 



# 

■ ■ 


■ 


■ 


8 
■ 


■ ■ 




I I 








c 


-z'/T 





s/cte view, converter cross-section 



a) side view, converter cross-section 



2/L 



m - shell thermocouple 
« - monolith thermocouple 



m - shell thermocouple 
a - monolith thermocouole 



r/Rmw=1 



r/Rmw=0 




shell thermocouple 
locations 



b) front view - converter cross-section 
Figure 4. Thermocouple Locations, Cases 1 and 2 



minor axis shell 
thermocouple locations 



major axis shell 
thermocouple locations 




b) front view, case 3 converter cross-section 
Figure 5. Thermocouple Locations, Case 3 
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RESULTS AND DISCUSSION 

PRESSURE LOSS - Table 5 lists predicted versus mea- 
sured pressure loss for all six cases. Pressure loss is 
underpredicted in cases 1 and 2, but overpredicted in 
case 3. There are several possible reasons for discrep- 
ancies in pressure loss predictions, including 

• uncertainty in monolith channel diameter; 

• uncertainty in the A/F ratio, which is used (along with 
fuel-flow measurements) to calculate the inlet mass 
flow rate of the CFD cases; 

• discretization and numerical error (e.g., numerical 
diffusion); 

• turbulence modeling inaccuracies; 

• transient effects of flow pulsations. 

Table 5: Pressure Loss Results 





Pressure Loss 


(kPa) 


case 


measured 


computed 


% difference 


1M 


10.8 


9.35 


-14.4% 


1H 


15.81 


12.53 


-23.1% 


2M 


11.71 


11.17 


-4.7% 


2H 


20.89 


18.56 


-11.8% 


3M 


8.63 


11.43 


27.9% 


3H 


16.18 


18.5 


13.4% 



CONVERSION - Table 6 lists conversion efficiencies for 
CO and total hydrocarbons for all six cases. Hydrocar- 
bon efficiencies are fairly accurate in all cases; carbon 
monoxide efficiencies are underpredicted in all cases. A 
probable reason for the consistent CO underprediction is 
the lack of NO x chemistry in the computations. In the 
experiment, significant NO x reduction takes place, which 
supplies further oxygen to continue the CO oxidation 
reaction. 

It should be noted that the catalytic converters in this 
study are not designed for use with the test engine, nor is 
the engine controlled for optimum use with the convert- 
ers. Therefore complete conversion is not expected. In 
all cases, it is lack of excess oxygen which keeps the pol- 
lutant species from completely oxidizing. 

Table 6: Conversion Results 



case 


species 


conversion effi 
measured 


ciency (%) 
computed 


2M 


CO 


71 


55 




HC-totat 


71 


74 


2H 


CO 


47 


33 




HC-total 


71 


64 


3M 


CO 


84 


59 




HC>total 


77 


70 


3H 


CO 


65 


40 




HC-total 


76 


69 



TEMPERATURE: CASES 1M AND 1H - Figures 6 and 7 
show measured versus computed temperature profiles in 
cases 1 M and 1 H, respectively. Because there is no reac- 
tion in cases 1M and 1H, no increase is seen in gas tem- 
peratures as flow moves through the monolith. In both 
cases, shell temperatures are significantly cooler than 
gas temperatures, especially in the monolith region. On 
the inlet cone, steep axial temperature gradients are 
seen, decreasing shell temperatures by several hundred 
degrees over the length of the cone. A similar profile is 
seen on the outlet cone. 

The temperature profile entering the monolith is non- 
uniform due to flow and temperature patterns in the inlet 
cone. Computed temperature contours showing this pat- 
tern at the monolith entrance are shown in Figure 8. A 
hot core temperature, which has not been significantly 
affected by thermal boundary layers in the inlet pipe and 
cone, forms the high-temperature core flow in the center 
of the monolith. Fluid which has been cooled somewhat 
by thermal boundary layers in the inlet pipe and cone 
enters the monolith near the outer radius. At the inlet 
monolith face (zYL* = 0), the predicted temperature differ- 
ence between the centerline and outer radius is approxi- 
mately 50 degrees C in cases 1M and 70 degrees C in 
case 1H. Similar temperature differences are seen in the 
measurements at z7L=.125. 

The lower-temperature region near the monolith outer 
radius grows along the length of the monolith until it 
reaches the monolith exit. In the monolith, the layer is 
best described as a conduction layer, as it is maintained 
by radial conduction through the insulating mat. Although 
the conduction layer continues to grow along the length 
of the monolith, it is more significantly impacted by the 
inlet temperature profile than by radial conduction in 
cases 1M and 1H. In case 1H for example, the tempera- 
ture difference (T c) - T mw ) which is a maximum at the 
monolith outlet, has already reached 60% of its maximum 
value at the monolith inlet. In cases 1M and 1H then, 
monolith temperature profiles are affected more signifi- 
cantly by flow in the inlet cones than by radial diffusion. 

The influence of the conduction layer is apparent in the 
measured and computed shell temperature shown in Fig- 
ures 6 and 7. Shell temperatures decrease to a mini- 
mum at approximately zYL* = .75, or 3 A along the length 
of the monolith. The decrease is due to the decreasing 
radial heat flux from the monolith as the conduction layer 
grows and radial temperature gradients in the monolith 
and mat become less severe. At around zYL = .75, axial 
conduction from the outlet cone shell begins to have a 
significant impact on shell temperature, and the shell 
temperature rises again. 
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Figure 6. Temperature Profiles, Case 1M 
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Figure 7. Temperature Profiles, Case 1 H 
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Figure 8. Temperature Contours in the Inlet Cone of 
Case 1M (side view, cross-section) 

TEMPERATURE.CASES 2M AND 2H - Figures 9 and 
10 show measured versus computed temperature pro- 
files in cases 2M and 2H, respectively. An increase in 
measured gas temperature of 62 and 37 degrees is seen 
in the monolith of cases 2M and 2H, respectively. Similar 
temperature gains are seen in the computations, as 
shown in Figures 9 and 10. This gas temperature 
increase is due to the heat of reaction for the three reac- 
tions considered. Under the high-temperature conditions 
considered here, the reactions take place over the first 
several centimeters of the front monolith, between z7L'=0 
and zVL-0.2. Methane oxidation is an exception; how- 
ever, this reaction generates an insignificant amount of 
heat compared to the other reactions considered. 

Note that there is more heat generated in case 2M than 
in 2H, although overall temperature levels are greater in 
case 2H. This behavior is due to the engine control 
scheme, which runs slightly rich at higher loads. Rich A/F 
ratios decrease 0 2 levels in the exhaust, leading in turn 
to less complete conversion in the higher-load 2H case. 
The same trend is seen in cases 3M and 3H. 

Unlike in case 1, in case 2 there is a significant difference 
between the monolith solid and fluid temperatures 
between z7L = 0 and zVU = .5. These monolith tempera- 
ture profiles in case 2H are shown in Figure 11. Solid 
temperature reaches a maximum immediately at the inlet 
face of the monolith. There is however a radial tempera- 
ture distribution across the monolith solid front face. In 
case 2H for example, [T d - T mw ] is in excess of 70 
degrees C. A significant amount of heat is transferred 
from the monolith solid to the fluid in the channels 
between zVU = 0 and zVL = 0.5 In this region, gas tem- 
peratures rise while solid temperatures fall, until the two 
become essentially equal around z7L'=0.5. 

The differences between solid and fluid temperatures in 
the monolith do not change the basic nature of the con- 
duction layer at the monolith outer edge. It is still influ- 
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enced heavily by gas inlet temperature profile as well as 
radial diffusion in the monolith itself. This conduction 
layer is shown in Figure 12, a plot of radial temperature 
profiles in the monolith at two axial locations in case 2H. 
At zVL-,125, both solid and fluid temperatures decrease 
significantly near the wall, although the solid profile is 30 
degrees or so hotter than the fluid profile. At z7L ,= .875, 
the solid and fluid temperatures are nearly identical. The 
significance of the monolith conduction layer and the exo- 
thermic reactions combine in Figure 12 to create an inter- 
esting fluid temperature field: centerline fluid temperature 
increases in the axial direction, while fluid temperature at 
r/R mw =1 decreases in the axial direction. 

The monolith temperature profiles shown in Figure 11 
lead to mat temperature maximums at the mat leading 
edge (at z'/L-O), adjacent to the monolith. It is here 
where the mat will first be pushed past its temperature 
limit in very-high temperature durability tests. As has 
been shown in case 2, this maximum mat temperature 
will be impacted significantly by chemical reaction and 
temperature distribution at the monolith inlet. Though a 
uniform flow distribution will improve backpressure and 
catalyst coating longevity, significant decreases in mat 
and shell temperatures can be realized by reducing flow 
in the outermost channels, i.e. as r/R mw approaches 1 . 
In the limit, elimination of flow in these channels would 
create a stagnant region without reaction or associated 
heat generation. Past work which has focused on the 
adverse effects of flow maldistribution in converter mono- 
liths (12,13) does not consider the durability effects of 
reduced flow at the far outer radius of the monolith. 

Shell temperature profiles between cases 1M and 2M 
look very similar. To better appreciate the difference 
between them, recall that the inlet gas temperature of 
case 2M is 30 degrees C cooler than case 1M. So the 
nearly equivalent computed shell temperatures in these 
two cases demonstrate the heating effect of the conver- 
sion reactions on shell temperatures. The temperature 
increases associated with chemical reaction are impor- 
tant when performing structural analysis of a catalytic 
converter assembly. Both mat properties and shell ther- 
mal stresses are heavily dependent on temperature. A 
major benefit of this type of analysis is its accurate pre- 
diction of mat and shell temperatures under reacting con- 
ditions, which can be used as input boundary conditions 
for structural finite element analysis. 
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Figure 9. Temperature Profiles, Case 2M 
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Figure 10. Temperature Profiles, Case 2H 
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Figure 1 1. Monolith Solid and Fluid Temperature 
Contours in Case 2H (side view, cross- 
section) 
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Figure 12. Monolith Conduction Layer Development in 
Case 2H 

TEMPERATURE: CASES 3M AND 3H - Figures 13 and 
14 show predicted and measured temperature profiles in 
cases 3M and 3H, respectively. Note that the shell tem- 
perature profile is significantly different from previous 
cases due to the different insulation strategy used in 
design 3 (reference Figure 1). Rapid axial decrease in 
shell temperature along the inlet cones is followed by a 
sudden rise at the mid-bed region, where the shell is not 
protected by insulation. 

In all cases of steep axial temperature gradients on the 
shell, it is imperative to include 3-D cells representing the 
steel converter shell. Axial conduction through these ele- 
ments controls temperature profiles in these critical high- 
solid-gradient regions. This axial conduction effect is cru- 
cial to the accurate prediction of mat-edge temperatures 
and shell temperature gradients, both of which are key 
predictors of converter durability. 



In cases 3M and 3H temperatures at the monolith outer 
edge, and in the mat and shell, are significantly cooler at 
the cross-sectional major axis than the minor axis. Mea- 
sured and computed temperatures are lower along the 
major axis than the minor axis by approximately 50C in 
case 3M and 70C in case 3H. This trend is seen every- 
where except at the outlet cone shell, where computed 
major and minor axis shell temperatures approach each 
other. 

This trend in shell temperatures can be traced to temper- 
ature profiles in the monolith. Cross-sectional monolith 
temperature profiles are shown in Figure 15. The differ- 
ences between major and minor axis temperature levels 
is due in large part to flow effects in the inlet cone. As in 
cases 1 and 2, significant decreases in gas temperatures 
are seen in the thermal boundary layers of the inlet pipe 
and cone. These boundary layers lead to decreased 
temperatures entering the monolith near its outer edge, 
especially along the major axis where there is more 
severe separation, recirculation, and temperature degra- 
dation in the inlet cone. A secondary factor is the flange 
on the converter shell, where the top and bottom shells 
are welded together. The flange acts as a heat sink 
which conducts heat away from the monolith along the 
major axis. 
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Figure 13. Temperature Profiles, Case 3M 

The previously identified monolith conduction layer is 
interrupted in the mid-bed region (between the monoliths) 
in cases 3M and 3H. Temperature contours in this region 
are shown in Figure 16. The gas which contacts the shell 
in the mid-bed region is the same gas which has been 
cooled below the core temperature in the monolith con- 
duction layer. Although heat transfer coefficients 
between the gas and mid-bed shell are high, the temper- 
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ature of the gas contacting the shell is lower than the 
centerline gas temperature by 100 degrees C or more. 
The shell is therefore protected from the maximum tem- 
peratures in the core flow. Although the conduction layer 
is interrupted by the gap between the monoliths, it contin- 
ues in the second monolith without significant changes. 
The leading edge of the second mat is therefore also pro- 
tected by the conduction layer in the mid-bed region. 
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Figure 14. Temperature Profiles, Case 3H 

The conduction layer in the mid-bed region has signifi- 
cant impact on converter durability. Chen et al. (5) have 
proposed that more insulation leads to a more uniform 
radial temperature profile in the front monolith, leading in 
turn to reduced thermal gradients and thermal stresses in 
the rear section of the monolith. On the other hand, a 
more insulated design retains higher temperatures in the 
conduction layer, leading to higher shell temperatures in 
the mid-bed region. These higher shell temperatures can 
lead to excessive shell warpage, thus affecting the dura- 
bility of the shell assembly. A tradeoff must be made 
which provides enough insulation to protect the monolith 
from transient thermal gradients without overheating the 
mid-bed shell. 

Temperatures are significantly lower in the second mat 
due to the expanded conduction layer in the second 
monolith. Though the second mat will be less affected by 
high-temperature degradation, it will be more susceptible 
to low-temperature problems should they arise. Note for 
example in case 3M, where inlet gas approaches 730 
degrees C, the shell over the second brick along the 
minor axis is not yet to the expansion temperature of the 
mat. Incomplete mat expansion can lead to problems in 
mat erosion and substrate retention. Low-temperature 
concerns are increasingly important as manufacturers 
move toward thicker mats in converter designs. 



Fluid which contacts the outlet cone shell has been 
cooled in the monolith conduction layer. Here, as in the 
mid-bed region, the conduction layer (now better 
described as a thermal boundary layer) offers the shell 
some protection from the hotter core flow. However, the 
experiments show hotter temperatures on the outlet cone 
surface than do the predictions. Thermal radiation from 
the substrate rear face (not accounted for in the computa- 
tions) may play a role in increasing cone shell tempera- 
tures, especially along the major axis where view factors 
between the two surfaces are highest. 
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Figure 15.Temperaure Profiles on Monolith Cross- 
sections, Case 3H (front view) 
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Figure 16. Temperature Profile in the Midbed Region of 
Case 3H (side view, major axis cross-section) 

CONCLUSIONS 

Significant conclusions of the study are as follows: 

1. A complete computational methodology has been 
established for prediction of catalytic converter ther- 
mal behavior under steady, high-load situations. The 
modeling method has been validated in 6 different 
cases. 

2. A thermal conduction layer, formed by entrance 
effects in the converter inlet cone and enhanced by 
radial heat diffusion in the monolith, has a significant 
impact on temperature behavior in the mounting mat, 
shell, and monolith, 

3. In the case of a dual-monolith, oval cross-section cat- 
alytic converter, significant variations in thermal 
behavior are seen along the major and minor axes. 
The mid-bed region also has a significant impact on 
shell and mat temperatures. 

4. The results obtained in this study may be used to 
develop equations for temperature field behavior as a 
function of design parameters such as brick radius, 
mid-bed size, shell thickness, mat thickness, etc. 
Such equations will be useful as design guidelines to 
be used throughout the catalytic converter design 
process. These temperature results will also be used 
as input boundary conditions for structural finite ele- 
ment analysis to predict shell stress and deflection. 

Currently, the heat transfer model is used extensively to 
assess the effects of particular design changes, such as 
upstream bends, cone geometry changes, etc., on dura- 
bility performance. Future work will include the develop- 
ment of expressions noted in conclusion 4. Work is also 
underway to further enhance the computational model, 
especially, to include the effects of thermal radiation 
between the substrate faces and the shell cones. 
Though the computational method for this radiation cal- 
culation is in place and has been implemented with some 
success, the emissive properties of a converter substrate 
face are not yet well known. 
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NOMENCLATURE 

VARIABLES 

a s : monolith surface area per unit volume (m 2 /m 3 ) 

a pt : surface area of platinum per unit volume (m 2 /m 3 ) 

C p : specific heat (kg/kJ) 

Cfji mole fraction of species i in fluid (mol/mol) 

c s/ : mole fraction of species i in solid (mol/mol) 

d h : hydraulic diameter (m) 

D mf - : mass diffusivity of species i in exhaust (m 2 /s) 

f : laminar friction factor (-) 

h: heat transfer coefficient (W/m 2 C) 

AH|: heat of reaction (J/mol,) 

k/ : reaction rate constant for reaction i (mol/m pt 2 s) 

m mj : mass transfer coefficient for species I (kgj/m 2 s) 

P : pressure (kPa) 

r : radial coordinate (m) 

Rj : specific reaction rate (mol/m 2 Pt s) 

S : convective source term (W/m 2 ) 

T : temperature (degrees Celsius) 

V : velocity (m/s) 

z : axial coordinate (m) 

OL : monolith porosity factor (-) 

X x ,Xy,X z : thermal conductivity in the x,y,z directions 

(W/mK) 

m: viscosity (kg/m 2 s) 
p: density (kg/m 3 ) 

NON-DIMENSIONAL PARAMETERS 

Re = pVjd/,/^: Reynolds number 

Nu = hd/j/Xfli Nusselt number 

Sh = m m /dfJD mi : Sherwood Number 

z + = z/(d/,Re): non-dimensional axial length 

SUBSCRIPTS 

cl : centerline 
f: fluid 

mw : monolith wall 
s : solid 

z : axial direction 

i : reaction identifier: i=1 - CO+(.5)0 2 =C0 2 

i=2 - C 3 Hq+(4.5)0 2 =3C0 2 +3H 2 0 
i=3 - CH 4 +30 2 =C0 2 +4H 2 0 
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